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YAWING  MOTION  OF  5.'0  MK  41  PROJECTILE 
STUDIED  BY  MEANS  OF  YAW  SONDES 


by 

W.  R.  Haseltine 
Research  Department 


ABSTRACT.  Hie  5'JO  Mk  4l  projectile  was  fired  at  both 
high  and  low  nuzzle  velocities,  and  its  oscillatory 
yawing  motion  studied  by  a  telemetry  technique.  At  high 
subsonic  velocities  this  projectile  develops  a  steady 
precesBional  motion  of  amplitude  about  4  degrees.  On  a 
long  flight,  from  high  velocity  launch,  this  behavior 
shows  up  late  on  the  downward  trajectory  leg.  The 
instability  at  small  yaw  and  the  limit  cycle  were  found 
to  be  caused  by  a  large  and  highly  nonlinear  Magnus 
moment.  Aerodynamic  coefficients  were  derived. 
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THE  PROBLEM 


In  the  past  fifteen  years  it  has  become  evident  that  many  spin- 
stabilized  projectiles  and  rockets  show  symptoms  of  dynamic  instability. 
While  they  are  gyroscopically  stable,  and  the  yawing  motion  is  oscilla¬ 
tory,  the  two  modes  of  oscillation  are  not  always  damped.  Over  certain 
speed  ranges  either  the  fast  mode,  the  nutation,  or  the  slow  mode,  the 
precession,  will  grow  from  small  amplitudes  to  larger  instead  of  decay¬ 
ing  away.  Some  otherwise  promising  designs  have  proved  worthless  for 
their  intended  purpose  because  of  such  behavior.  Even  among  those 
missiles  which  acre  well  enough  behaved  to  have  been  adopted  for  service 
use,  many,  perhaps  most,  show  an  instability  of  the  precession  at  high 
subsonic  speeds,  say  between  Mach  numbers  0.7  and  1.0.  Typically,  as 
in  the  case  of  the  105-mm  Shell,  Ml,  the  precession  grows  from  a  small 
size  to  an  amplitude  of  a  few  degrees,  and  then  stays  in  a  limit  cycle 
(Ref.  l).  Worse  behavior  is  also  known. 

A  few  years  ago,  when  development  of  rocket-assisted  ammunition  for 
the  Navy's  5"0  guns  was  undertaken,  it  became  important  to  know  whether 
the  5'JO  Shell,  Mk  4l  (that  for  the  5"0/54  gun)  showed  any  such  behavior; 
first,  because  the  very  long  flight  times  of  RAP  and  the  expected  low 
ratio  of  spin  to  velocity  late  in  flight  should  accentuate  the  pathology 
and  its  effects  on  range  and  dispersion:  second,  because  existing  spark 
range  data  (Ref.  2),  though  showing  considerable  scatter,  indicated  in¬ 
stability  at  high  subsonic  speeds.  New  spark  range  measurements  were 
made  by  the  Ballistic  Research  Laboratories  and  wind-tunnel  measurements 
by  the  Naval  Ordnance  Laboratory,  White  Oak.  Thes'~  measurements  and 
their  significance  are  reviewed  in  Ref.  5.  They  confirm  the  indications 
in  Ref.  2  that  the  Magnus  moment  coefficient  at  high  subsonic  speeds  and 
very  small  yaw  angles  is  such  that  undamping  of  the  precession  is  to  be 
expected.  In  addition,  the  wind-tunnel  data  indicate  a  change  of  sign 
of  the  Magnus  moment  at  about  4  degrees  yaw.  L'uch  a  change  points  to 
the  existence  of  a  limit  cycle. 

It  is  plainly  of  great  interest  and  importance  to  confirm  the  pre¬ 
dicted  behavior  by  experiments  on  real  shell,  observing  the  yawing 
behavior  over  whole  trajectories.  This  is  especially  true  since 
(a)  the  spark  range  experiments  were  limited  in  the  number  of  speed  and 
yaw  amplitude  combinations  obtained,  the  values  of  the  Magnus  moment 
and  the  yaw  damping  coefficients  scattered  badly  in  the  realm  of  inter¬ 
est,  and  there  were  insufficient  data  to  even  begin  to  deduce  the  func¬ 
tional  form  of  the  dependence  on  yaw  angle;  (b)  wind-tunnel  measurements 
of  Magnus  moment  at  small  angles  are  extremely  difficult  and  at  the  time 
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of  Ref.  3  did  not  exist  for  this  round  below  4  degrees,  and  there  are 
serious  questions  concerning  the  effects  of  Reynold's  number  and  the 
supporting  sting.  What  is  needed  is  a  device  which  will  do  for  shell 
what  the  Solar  Yaw  Camera  (Ref.  4,  p.  405)  did  for  rockets,  and,  if 
possible,  easier. 

Now  even  before  the  problem  became  acute  for  the  5'*0/54,  the 
British  had  announced  (at  the  second  meeting  of  Panel  07,  TTCP,  in  1963) 
that  they  had  developed  a  Yaw  Sonde  with  Just  the  desired  properties. 
Their  government  made  a  gift  of  a  few  sondes  to  this  Station  in  1964, 
and  we  later  purchased  more.  Use  of  the  sonde  by  its  originators  is 
described  in  Ref.  5-  It  was  therefore  decided  to  use  these  sondes  to 
study  the  motion  of  the  Mk  4l  shell.  Two  groups  of  shots  were  planned. 
One  group  was  to  be  fired  with  a  muzzle  velocity  just  above  sonic,  so 
that  a  long  flight  at  speed  between  Mach  No.  0.7  and  1.0  would  be  ob¬ 
tained,  giving  full  opportunity  to  observe  limit  cycle  behavior.  The 
other  group  was  to  be  fired  at  full  service  charge,  and  for  maximum 
range,  to  check  the  theoretical  prediction  that  the  round  would  fly  at 
small  yaw  until  several  seconds  before  impact,  when  a  limit  cycle  in 
precession  would  develop.  Later,  these  plans  were  extended  to  include 
tests  of  rounds  without  the  rotating  band,  for  the  results  reported  in 
Ref.  6  suggested  that  decided  changes  might  result.  All  records  of 
yawing  motion  were  to  be  fitted  with  theoretical  curves  calculated  by  a 
computer,  so  that  values  could  be  derived  for  the  round's  aerodynamic 
coefficients,  including  insofar  as  possible  their  functional  dependence 
on  yaw  angle  (see  Ref.  7).  It  was  hoped  that  this  set  of  experiments 
would  be  the  start  of  a  series  studying  more  generally  Magnus  instabili¬ 
ties  of  spinning  projectiles  and  trying  to  find  cures. 


PRINCIPAL  CONCLUSIONS 


The  5-0  Shell  Mk  4l  does  indeed  exhibit  a  limit  cycle  in  precession 
at  Mach  numbers  between  0.7  and  about  1.0.  When  it  is  fired  at  a  veloc¬ 
ity  just  above  sonic,  the  small  initial  precession  starts  to  increase 
almost  at  once,  and  then  reaches  a  steady  amplitude.  In  all  but  two  of 
the  nine  cases  available  this  amplitude  was  a  little  less  than  4.0 
degrees,  and  there  is  some  doubt  as  to  the  significance  of  the  two  excep¬ 
tions.  The  initial  nutation  is  damped,  at  least  at  first,  but  may  per¬ 
sist  at  small  (1.0  deg  or  less)  amplitude  superimposed  on  the  larger 
precession.  When  the  shell  is  fired  at  about  Mach  2.2  and  a  Q.E.  of 
45  deg.,  both  precession  and  nutation  are  rapidly  damped  out  to  undetect¬ 
able  size.  The  flight  time  is  about  74  sec,  and  at  from  15  to  25  sec 
before  impact  small  precessions  become  visible,  which  then  grow  steadily, 
but  do  not  reach  the  4.0-deg  amplitude  before  impact.  Nutations  are 
invisible  or  nearly  so  during  this  time. 
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In  the  comparison  of  computed  curves  of  orientation  angle  with  the 
data  from  the  sonde s,  very  good  fits  were  obtained  to  both  frequency  and 
amplitude  of  the  precessions,  and,  when  the  data  were  of  adequate  quality 
to  show  nutations  reliably,  also  to  the  frequency  of  these.  It  appears 
certain  that  from  Mach  number  0.70  to  at  least  0.95  and  possibly  1.00, 
the  Magnus  moment  coefficient  ICp  (or  -ir/l6  C  )  is  positive  and  at 
least  as  large  as  0.095,  quite  possibly  as  largeas  0.150;  and  that  it 
changes  sign  rather  suddenly  at  4.0  degrees  or  a  little  more  total  yaw 
angle,  staying  negative  at  least  to  6.0. deg.  Supersonically  the  value 
of  was  assumed  to  be  -0.024,  and  the  fits  are  all  relatively  insensi- 

tiveito  moderate  changes  in  this  value.  The  yaw  damping  factor 

kh  ( [c  +  C^]  )  was  in  most  of  the  comparisons  taken  to  have  a 

fixed  value,  since  the  nutations  were,  except  for  one  case,  too  poorly 
defined  for  amplitude  fitting.  The  value  chosen  was  5-0.  It  is  this 
value  which  together  with  a  K^,  of  0.095  gave  a  good  fit  to  precessional 
growth  rate  in  the  subsonic  firings.  The  absence  of  noticeable  persist¬ 
ent  nutations  in  the  computed  curves  strongly  suggests  that  this  value 
of  K„  is  too  large,  and  that  therefore  the  K_  is  too  small.  Results 
for  tne  one  round  with  an  extended,  clean  record  of  nutations  appear  to 
indicate  that  KH  Is  quite  small  between  Mach  0.9  and  1.0,  but  larger 
below  0.9-  There  is,  however,  no  indication  that  Kg  is  ever  negative. 

In  order  to  secure  good  fits  in  frequency,  it  was  necessary  to 
assume  that  the  main  moment  coefficient  (tt/8  C^)  was  a  function  of 

both  Mach  number  and  the  square  of  the  yaw. 

No  evident  difference  was  found  between  the  behaviors  of  the  three 
varieties  of  the  shell  tested.  Neither  machining  the  rotating  band 
flush  with  the  rear  bourrelet  nor  removing  the  band  entirely  seemed  to 
make  an  appreciable  difference  in  flight  behavior,  provided  always  that 
muzzle  velocity  and  Initial  spin  were  held  the  same. 

A  consistent  peculiarity  in  the  data  from  the  supersonic  shots  is 
attributed  to  a  strong  wind  shear  between  9,000  and  11,000  ft.  altitude. 


THE  SONDE 

The  "yaw  sonde"  is  essentially  a  miniature  radio  transmitter  which 
is  frequency  modulated  by  a  photocell,  this  photocell  being  periodically 
illuminated  by  the  sun  shining  through  a  pinhole.  The  pinhole  Is  in  the 
periphery  of  the  round.  The  photocell  is  a  circular  wafer  of  silicon 
mounted  parallel  to  the  round's  axis  and  located  between  the  pinhole  and 
axis.  The  face  of  the  disk  is  masked  except  for  two  narrow  strips 
arranged  in  a  "V"  pattern.  In  flight,  once  each  revolution  a  small 
image  of  the  sun,  formed  by  the  pinhole,  sweeps  across  the  disk,  hitting 
first  one  arm  of  the  "V",  then  the  other.  The  position  along  the  lengths 
of  the  slits  at  which  it  crosses  depends  on  the  angle  between  the  axis 
of  the  round  and  a  line  to  the  sun. 


5 


In  the  original  version  of  the  sonde,  the  photocell  disk  was  mounted 
on  the  axis,  and  the  field  of  view  was  about  20°.  In  order  to  increase 
this  acceptance  angle,  so  as  to  permit  more  flexibility  in  firing  con¬ 
ditions,  this  Station  moved  the  cell  closer  to  the  pinhole,  filled  the 
space  between  pinhole  and  disk  with  clear  plastic,  and  added  a  transis¬ 
tor  preamplifier.  An  acceptance  angle  of  a  little  more  than  90°  was 
obtained,  at  some  sacrifice  in  resolution.  In  Fig.  1  are  shown  the 
transmitter,  which  is  contained  in  a  fuze  body,  and  the  packaged  cell 
and  preamplifier.  The  latter  is,  in  use  attached  to  the  inner  side  of 
the  cylindrical  portion  of  the  shell,  just  under  the  pinhole.  The 
battery  used  was  a  standard  fuze  type  reserve  battery,  housed  in  the 
rear  part  of  the  fuze-shaped  body. 

The  radio  frequency  signal  from  the  sonde  is  received  by  a  ground 
station,  and  the  demodulated  version  recorded  on  magnetic  tape  along 
with  timing  pulses  In  I.R.I.G.  format  B.  This  demodulated  signal  con¬ 
sists  in  a  train  of  pulse  pairs,  one  pair  for  each  revolution.  Oscillo¬ 
scope  traces  of  sections  of  two  such  records  are  shown  in  Fig.  2,  the 
upper  one  a  good  and  useful  signal,  the  lower  at  or  below  the  margin  of 
uselessness.  The  time  between  centers  of  pulse  pairs  is  the  reciprocal 
of  the  spin  frequency,  and  the  ratio  of  the  time  between  pulses  of  a 
pair  to  that  between  pairs  is  a  unique  measure  of  the  angle  sought — that 
between  the  round  axis  and  a  line  to  the  sun.  After  each  sonde  was 
mounted  In  the  projectile,  but  before  potting,  careful  measurements  were 
made  of  the  functional  relation  between  this  ratio  and  that  angle;  and 
this  relationship  was  later  used  in  reducing  the  data. 

The  data  on  the  tape  were  converted  to  digital  data,  also  on  tape, 
by  NODAC,  under  human  supervision.  This  new  tape  in  turn  was  subjected 
to  processing  by  two  successive  digital  computer  programs,  the  final 
result  being  a  tabulation  of  spin  rate  and  orientation  angle  and  its 
sine  against  time.  Plots  of  the  data,  usually  even  more  useful  than  the 
tables,  were  also  made. 


THE  ROUND 


All  the  experimental  projectiles  used  started  out  as  empty  standard 
5'.'0  Mk  4l  shell.  On  six  of  them  the  rotating  band  was  turned  flush  with 
the  rear  bourrelet,  and  on  four  the  band  was  split  and  entirely  removed. 
Each  projectile  was  then  sawn  In  half,  transversely.  The  front  portion 
was  turned  internally  to  a  smooth  cylindrical  surface  of  length  about 
375,  parallel  to  the  axis,  and  then  threaded  for  0775.  At  about  2725 
forward  of  the  cut,  a  conical  hole  (90°  full  angle)  was  machined,  the 
axis  of  the  cone  being  radial  with  respect  to  the  cylindrical  shell 
body.  The  intersection  of  the  tip  of  this  cone  with  the  inner  machined 
cylindrical  surface  formed  the  pinhole,  about  07050  diameter.  Ahead 
and  behind  this  depression,  two  countersunk  holes  were  drilled.  The 
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I  K'..  I  British  Yaw  Sonde,  With  Sensor  Modified  by  Naval  Weapons  Center. 
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FIG.  2.  Oscilloscope  Traces  of  Yaw  Sonde  Signals. 
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after  part  of  the  projectile  was  also  threaded,  from  front  to  back,  and 
an  externally  threaded  sleeve  was  used  to  connect  the  two  parts  back  to¬ 
gether  (eventually).  A  non-threaded  spacer  ring,  slipping  tightly  over 
the  threaded  sleeve  was  used  to  bring  the  reassembled  projectile  back  to 
proper  length.  Because  the  sonde  transmitter  and  battery  case  were 
housed  in  a  body  conforming  to  the  standard  NATO  artillery  fuze  (Army 
M51  family  shape),  and  threads,  an  adapter  was  required  between  the  shell 
body  and  sonde  body.  This  provided  proper  threads,  length,  and  contour. 

While  the  shell  body  was  still  in  two  pieces,  the  assembly  of  photo¬ 
cell  and  preamplifier  was  attached  by  screws  to  the  interior  of  the 
front  portion,  just  under  the  pinhole.  Then  the  sonde  body  with  trans¬ 
mitter  and  battery  (with  adapter)  was  screwed  in  and  electrical  connec¬ 
tions  made.  The  sonde  was  calibrated,  using  an  external  source  of  volt¬ 
age,  a  meter,  the  sun  (vie  a  heliostat)  as  a  light  source  and  a  dividing 
head.  After  the  rear  and  front  parts  had  been  reassembled,  the  interior 
was  filled  through  the  base  orifice  with  an  epoxy  potting  compound,  and 
the  base  plug  inserted.  Figure  5  is  a  photograph  of  two  rounds,  one  with 
the  standard  rotating  band  and  one  with  none. 

It  is  obvious  that  the  rounds  with  modified  or  no  rotating  band  had 
to  be  provided  with  sabots  to  provide  obduration  and  to  impart  spin. 

The  sabot  designed,  Fig.  k,  had  an  aluminum  body,  a  copper  skirt,  (nearly 
split),  and  two  steel  pins  which  fitted  into  holes  drilled  in  the  projec¬ 
tile  base.  This  was  very  successful. 

A  special  experiment  was  made  to  determine  the  powder  charges  needed 
to  reacn  the  desired  muzzle  velocities  (1240  and  2^50  ft/sec)  with  both 
standard  and  saboted  rounds.  The  2^50  ft/sec  was  chosen  deliberately 
less  than  service  M.V.  to  ensure  safety  with  saboted  rounds. 

Table  1  (Appendix  A)  contains  a  list  of  the  physical  properties  of 
all  rounds  for  which  aerodynamic  analysis  was  performed. 


hANGE  INSTRUMENTATION  AND  THE  FIRINGS 


The  signals  from  the  yaw  sondes  were  received  by  two  antenna  arrays, 
one  vertically,  the  other  horizontally  polarized.  Each  fed  a  standard 
automatic  frequency-tracking  receiver  and  a  pan  receiver  with  scope.  The 
latter  permitted  rapid  location  of  the  signal  and  lock-or. .  The  demodu¬ 
lated  output  of  each  receiver  was  recorded  on  tape,  both  direct  and  via 
frequency  modulation.  IRIG  timing  signals  were  recorded  on  the  same  tape. 

To  obtain  trajectory  and  velocity  history,  the  projectiles  were 
tracked  with  an  MSG3A  radar,  modified  to  provide  improved  tracking  in 
range  and  digital  output  data. 
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FIG.  3  5*0  Mk  41  Shell  With  Yaw  Sondes 
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Mk  41  Shell  With  Modified  or  No  Rotating  BanJ. 


f 

i 


NVC  TP  4779 


A  modified  Hawk  radar  was  used  in  an  attempt  to  get  better  veloci¬ 
ties  by  the  Doppler  method.  Some  good  signals  were  obtained  but  they 
always  started  too  late  to  be  of  much  use. 

The  initial  velocities  were  measured  with  sky  screens,  main  reliance 
being  placed  on  a  pair  built  locally  and  located  at  150  and  200  ft  down- 
range  on  the  gun  line.  They  operated  almost  perfectly.  Backup  was  pro¬ 
vided  by  a  French  2-meter-base  sky  screen  on  loan  from  Aberdeen  Proving 
Ground,  and  by  Bowen  cameras.  The  axes  of  the  sky  screens  were  aligned 
perpendicular  to  the  gun  bore. 

Standard  meteorological  data  were  obtained,  surface  temperature 
pressure  and  wind  being  made  at  about  quarter  hour  intervals  during  each 
test,  and  a  balloon  ascension  to  obtain  corresponding  upper-air  data  was 
made  shortly  before  commencement  of  each  day's  firings. 

Prior  to  each  day's  shooting,  the  time  of  start  and  the  azimuth  of 
fire  were  carefully  chosen  to  maximize  the  time  that  the  sun  was  in 
position  to  be  within  the  receptive  cone  of  a  sonde  during  a  whole 
trajectory. 

The  first  test,  E-6293  on  15  December  1966,  showed  that  the  methods 
devised  worked,  end  that  one  could  hope  to  get  good  yaw  records  from  half 
or  more  of  the  rounds  fired.  Radar  coverage  was  marginal.  It  was  hoped 
to  complete  the  program  of  firing  by  the  end  of  calendar  year  1967. 
Crowded  schedules,  bad  weather,  and  repeated  severe  difficulties  with 
radar,  sky-screens  and  timing  led  to  delays  and  a  series  of  abortive 
tests,  so  that  the  shoots  weren't  completed  until  14  May  1968.  Diffi¬ 
culties  with  data  reduction  and  illness  of  the  author  delayed  the  analy¬ 
sis  of  the  data  by  curve  fitting  until  early  1969. 

Table  2  is  a  list  of  all  rounds  for  which  yaw  sonde  signals  were 
received,  with  firing  condition.  Table  3  lists  pertinent  meteorological 
data. 


THE  THEORY 


The  equations  of  motion  used  in  the  computation  of  the  curves  which 
were  fitted  to  the  data  may  be  divided  into  two  sets.  The  first  set 
describes  a  planar  particle  trajectory  ar.d  a  simple  spin  history.  It 
also  includes  an  equation  for  generation  of  an  exponential  air  density. 
The  second  set  describes  the  angular  motion  of  the  projectile.  It  is 
linearized  in  that  all  angles  both  of  orientation  and  direction  of 
motion  are  treated  as  small  deviations  from  those  appropi 'ate  to  a  pro¬ 
jectile  following  with  perfect  trail  a  particle  trajectory.  This  still 
permits  nonlinearities  to  appear  in  the  aerodynamic  force  system.  Also 
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included  in  the  second  set  is  a  calculation  of  the  angle  represented  by 
the  sonde  data.  A  list  of  symbols  is  provided  in  Nomenclature  section. 

The  force  and  moment  systems  used  may  be  described  (except  for  drag) 
in  the  body  fixed  axis  system  of  Ref.  8.  The  formulas  are  given  here  in 
both  aerodynamic  and  ballistic  notation.  The  first  has  the  advantage  of 
familiarity  to  a  wide  circle,  but  the  disadvantage  of  some  ambiguity  or 
lack  of  agreement  concerning  sign  conventions  and  use  of  some  factors  of 
2  in  the  definition  of  the  derivative  coefficients.  The  second  is  less 
known,  but  is  completely  unambiguous. 


FY  =  QSCaa  v/v  =  "  p  d\Vv 

F,  =  QSCza  w/V  -  -  p  d\vw 

M*  =  QSd  (pd/2V)  =  -  pd4KA  Vp 

"y  -  [C«  +  Cnpe  (pd/2V)v/»4(0m<i*Cn>i)(qd/2V)] 

■  Dd5  [k^Vv  -  Hj  pvd  -  KjjVqdJ 

MZ  *  Ka  <»/*)  +  C„p6  ( pd/2V ) v/V  +  ( 

+  Cna)(rd/2V) J 

=  pd'  [  -  l^Vv  +  K^,  pvd  -  K^Vr  d  J  . 

The  drag  force  is  taken  as  directed  along  the  relative  wind: 
PD  .  QC  CD  =  °  d2  V2  Kd  . 

This  set  of  definitions  implies  the  following  relations: 


kd  .  (ir/3)  CD 

Ks  ■  '  (w/8>  cm 

ka  =  -  (ir/l6) 

"h  ■  W0> 
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K,,  -  -  (ir/16)  Cnp? 

h  ’  '  *  <W 


In  principle  all  the  coefficients  are  functions  of  Mach  number,  V/a, 
and  of  the  square  of  the  yaw,  |b|  .  In  practice  only  certain  ones  were 
so  taken: 


Kj  -  fen  (I&I2)  *(m)  +  Ktpg  (1-f)] 


In  both  of  these  forms  |&|  has  its  computed  instantaneous  value.  The 
function  f(n)  was  almost  arbitrarily  taken  to  be  equal  to  unity  for 
U  <  1.0,  to  decrease  linearly  in  p  to  zero  at  )i  <  1.1,  and  to  stay 
zero  for  4  >  1.1. 

h  -  ho  *  KI»2  •  N5  • 


Here  j& is  a  smooth  function  of  the  time,  chosen  in  advance  of  each 
particular  computation,  and  intended  to  represent  a  sort  of  average  of 
the  running  |b|  .  K^g  itself  was  obtained  as  k^  K^qq,  with  KDQ(u) 

shown  graphically  in  Fig.  5,  for  which  an  abbreviated  listing  is  provided 
as  Table  4.  is  very  nearly  the  drag  coefficient  of  Refs.  2  and  9. 

KyQ  was  inserted  as  kg  Kj^fu)  where  is  as  follows: 


kM00 


u  <  0.755 
0.755  <  u  <  0.910 
0.910  <  u  <  0.96 
0.96  <  u  <  1.10 

p  >  1.10 

This  form  was  chosen  as  a  result 


1.621 

increases  linearly  to  1.957 
decreases  linearly  to  1.752 
constant  at  1.752 

constant  at  1.752  for  subsonic  shots; 
decreases  linearly  for  supersonic 
shots,  being  1.522  at  p  »  2.00. 

of  several  preliminary  fitting  attempts. 
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For  convenience,  certain  auxiliary  quantities  are  defined: 
U  -  V/d 
v  *  2rrn  /v 
D  =  pd3  Kjj/b 
N  =  pd3  Kjj/m 
L  *  N  -  D 
A  s  Pd5  KA/Ip 

M  =  pd5 

H  =  pd5^/^ 

T  =  p  d5  l^/lp 
G1  *  g  sin  a/dU2 
G2  =  g  cos  a/dU2  . 


The  independent  variable  corresponding  to  machine  time  was  chosen 
to  be  8,  the  dimensionless  arc  length.  This  variable  never  appears 
explicitly  in  the  equations,  serving  only  as  a  variable  of  integration. 
Its  time  derivative  is  U. 

The  first  set  of  equations  is 


t' 

=  l/u 

(1) 

U' 

«  -  DU  -  g  sin  a/Ud 

(2) 

a' 

g 

•  -  Gg  «  -  g  cos  a/d  U 

(3) 

n* 

*  -An 

(*) 

p' 

.  -  ydp 

(5) 
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x'  * 

d 

cos 

a 

(6) 

yl  - 

d 

sin 

a 

(7) 

The  second  i 

set 

of  equations  is 

6  = 

i 

-  e 

(8) 

0  = 

L 

&  + 

Gl9 

(9) 

0"  = 

[i 

L  V- 

(H  -  D  -  G.  )]  0'  +  (M  -  i  VT)  6 

X 

(10) 

+  1VG2 


sin  a  =  [(cos  a  -  0_  sin  a)  cos  A^  +  0,  sin  A.,]  cos  E 

2  S  3  S  S  (n) 

+  (sin  a  +  02  cos  a)  sin  Eg 

Equations  8,  9,  and.  10  each  represent,  of  course,  two  real 
equations.  Note  that  "sin  cr"  is  really  the  cosine  of  the  angle  between 
the  projectile  axis  and  a  line  to  the  sun. 


Note  that  no  provision  was  made  to  take  account  of  wind. 


The  equations  were  realized  on  an  Electronics  Associated  computer. 
Almost  all  components  available  from  one  console  and  about  half  those 
of  another  were  used. 


K,p2  was  taken  as  constant,  approximately  -.024.  KTL( | 6 | 2 )  was  set 

in  as  a  series  of  three  or  four  segments  each  linear  in  |&|2.  The 
general  form  chosen  was  based  on  previous  experience  with  other  projec¬ 
tiles,  and  proved  to  be  quite  adequate  in  this  case. 


CURVE  FITTING 


Samples  of  the  data  and  fits  to  them  are  shown  in  Appendix  B,  Plates 
IA-VIE.  For  example,  Plate  IA  is  a  plot  of  velocity,  trajectory  and  spin 
data  for  Round  11  of  E-6295,  and  the  overlay  IB  shows  the  corresponding 
computed  curve.  Plate  IC  is  the  yaw  sonde  record,  as  sin  <r,  while  the 
overlay  ID  contains  the  computed  version,  as  well  as  |b|  ,  pQT  (  |6p)/p 

and  Taking  the  vertical  distance  between  a  *  0.0  and 

a  *  1.0  as  unity  for  the  other  curves  also,  the  plots  represent  100  |&|2 

in  radians  squared,  1.2  *  10^T,  and  0.457  In  this  particular  case 

PQ  T/p  is  5  •  10  5  for  &2  small,  and  is  zero  at  62  *  50*0  “  10  ^  rad2. 
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The  smooth  curve  approximating  |6 is  the  |6 |^(t)  that  was  used  in 
this  case.  The  other  sets  of  plates  correspond,  in  general,  to  the  same 
pattern . 


Fitting  for  the  subsonic  shots  proceeded  as  follows.  A  plot  of  V, 
y,  x  and  n  versus  time  was  placed  on  the  table  and  overlaid  with  vellum. 
The  physical  parameters  of  the  round,  surface  air  density,  velocity  of 
sound,  and  guessed  values  of  the  aerodynamic  coefficient  functions  were 
inserted  (as  k.,  k2,  T  (62),  N,  etc.)  by  potentiometer  settings,  as  were 
also  initial  conditions  U  ,  a  ,  n  ,  and  computed  curves  of  the  four 
dependent  variables  drawn  by  tne  plotter,  driven  by  the  computer.  By 
successive  trials  of  different  k.  (drag)  values  and  of  slight  changes 
in  a  and  occasionally  U  ,  the^fits  between  computed  and  measured 
curves  of  V  and  y  were  improved.  Then  the  n  fit  was  improved  by 
adjusting  Ka» 


The  second  step  (after  setting  in  sun  azimuth  and  elevation)  was  to 
place  a  plot  of  the  measured  sin  a  vs.  t  curve  on  the  plot  table, 
overlay  it  with  vellum,  guess  initial  0~,  01  and  0i,  (©2  and  ©, 
were  taken  always  as  zero  initially),  ana'have  the  machine  superimpose 
its  computed  sin  a  curve  on  the  data.  Then  the  sun  angles,  Ag  and  Eg 
were  adjusted  by  small  amounts  and  by  tried  and  error  until  the  general 
trends  of  the  two  sin  a  curves  were  in  fair  agreement.  This  adjustment 
was  usually  1  to  3  degrees,  but  in  two  cases,  that  to  Eg  was  much  larger, 
and  in  one  case  that  to  Ac.  Following  this,  adjustments  were  made  to  the 

T(6  )  function,  and  occasionally  to  kg  (affecting  Kj^).  Emphasis  at  this 
stage  was  on  ensuring  that  there  were  no  major  discrepancies  in  the  com¬ 
puted  yawing  frequencies,  and  that  the  initial  amplitude,  rate  of  growth 
and  final  steady  amplitude  of  the  precessions  were  nearly  right. 


Since  it  had  been  found  in  early  trials  that  good  fits  to  both  the 
V  and  y  curves  could  not  be  obtained  without  allowing  some  effect  of 
the  yaw  on  drag,  the  computed  62  was  plotted,  and  a  function  unit  set 
up  tc  provide  a  smooth  function  of  time  nearly  reproducing  a  running  mean 
of  62,  Attention  was  then  returned  to  the  plot  of  trajectory  variables, 
k^  and  aQ  were  readjusted,  now  along  with  the  parameter  until 

satisfactory  matches  in  V  and  y  were  obtained  simultaneously,  (when 
this  had  succeeded,  the  x  curves  were  also  satisfactory.)  Any  necessary 
readjustment  of  was  then  made. 

Finally  a  definitive  fit  to  the  sin  a  curve  was  made,  again  by 
trial -and-error.  This  time  K  was  varied,  as  well  as  k0  and 

P  IDO 

T(6  ),  and  the  initial  conditions.  Here  an  attempt  was  made  not  only  to 
reproduce  the  amplitude  behavior  of  the  precessions  but  also  their  phase. 
We  tried  to  obtain  good  fits  to  phase  and  frequency  both  on  the  very 
early  and  very  late  oscillations  and  a  good  average  frequency  in  the 
middle.  No  phase  errors  over  about  l/k  cycle  are  evident  in  the  final 
fits. 


I 

I 

I 

1 

I 

I 

3 
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Control  of  Magnus  moment  was  exercised  through  the  function 

T0(62)  .o0T/b  -P0i5Hr/lp  • 

2  2  2 

For  most  fits  this  was  a  constant  from  6  *  0  to  6  «  6.  ,  typically 

-4  2  ^ 

40.0  •  10  rad  ,  then  decreasing  linearly  to  some  negative  value  at 

2  2  2  2  -4 

6  *  6p,  becoming  zero  at  8  *  6,  +  10.0  •  10  ,  and  being  a  negative 

2  2  1 

constant  for  6  >  6^  .  In  a  few  cases  this  functional  behavior  was 

replaced  by  one  breaking  more  gently,  i.e.,  positive  TQ  »  Tqq  >  0  for 

62  •  10^  <  37.0,  decreasing  linearly  to  0.80  T  at  &2  .  10^  *  42.0, 

2  00  Ok 

then  linearly  to  a  negative  value  at  8p  ,  being  zero  at  5  •  10  ■  50-0, 

and  a  negative  constant  for  6  >  &2  .  The  substitution  of  this  function 

for  the  previous  one  (for  the  &2  •  10^  *  40.0)  made  no  change  which  could 

be  detected  in  the  yawing  motion.  Broadening  the  transition  further,  by 

2 

making  the  decrease  from  Tqq  to  zero  take  place  linearly  over  a  6^ 

range  of  20  •  10  was  less  satisfactory.  Increasing  the  steepness  of 
the  transition  helped  the  behavior  of  nutation  amplitude  slightly,  but 
could  not  be  tolerated  during  much  of  the  time  for  it  caused  erratic 
computation.  Most  of  this  erratic  behavior  was  eventually  traced  to  a 
particularly  touchy  electronic  multiplier,  but  too  late  for  new  trials 
to  be  attempted.  For  all  rounds  but  numbers  5  and  7  of  E-7274, 

3  April  1968,  it  was  satisfactory  to  place  the  zero  of  T  at 

2  — li  2 

5  *  50.0  •  10”4  rad  .  For  the  exceptional  two  rounds  it  was  necessary 

2  2 
to  shift  the  whole  T(&  )  curve  bodily  along  the  6  axis. 

In  a  few  cases  the  early  nutations  were  clear  enough  so  that  it  was 

worth  attempting  to  fit  their  phase  and  frequency.  This  was  tried  using 

plots  with  an  enlarged  time  scale.  Small  adjustments  of  l£  or  less  were 
needed  to  the  nominal  ratio  I  /i^.  ,  and  corresponding  changes  in  4g(l^) 

to  secure  fit  of  both  nutation  and  precession  frequencies.  These  adjust¬ 
ments  are  well  within  the  precision  of  measurement  of  that  ratio.  The 
fit  so  obtained  was  the  best  one  for  the  first  few  seconds.  On  returning 

to  the  plot  for  the  whole  flight,  using  the  new  I  /i  ,  a  further  small 

P  * 

adjustment  of  was  needed  to  get  the  best  phase  behavior  averaged  over 

the  whole  trajectory.  In  one  case  the  nutations  were  clear  enough  so 
that  comparison  of  computed  with  measured  amplitudes  was  meaningful. 

In  this  case  it  was  not  possible  to  fit  these  amplitudes  and  the  pre¬ 
cession  behavior  using  a  constant  In  fact,  in  order  to  get  the  nuta¬ 

tions  to  persist  long  enough  for  a  10-second  fit  in  phase,  a  combination 
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of  and 

in  the  large  was  considered.  Implications  are  discussed  in  a  later 
section. 

Kjg  was  taken  as  the  constant  0.80  for  all  rounds.  It  was  intended 

that  Kjj  would  be  5.00  for  all  rounds.  By  inadvertence  some  of  the 

actual  values  deviated  a  little  from  this,  but  by  amounts  which  (as  shown 
by  actual  test  in  the  worst  case)  did  not  visibly  affect  the  plots. 

The  data  from  the  supersonic  rounds  were  fitted  in  the  same  way, 
except  for  the  four  following  changes: 

1.  Kp  was  not  made  to  depend  on  |b2|  .  Instead,  a  constant  was 
added  to 

2.  The  function  K^u)  was  different  for  |i^  1.10.  (This  change 

would  not  have  .affected  the  subsonic  fits  except  to  require  slightly 
different  initial  conditions,  for  on  none  of  these  was  p  >  1.1  where 
m  there  were  sonde  data.) 

3.  To  obtain  the  late  precessions  in  anything  like  the  measured 
amplitude,  and  without  using  combinations  of  H  and  T  which  were  very 
different  from  those  giving  fits  to  the  subronic  data,  it  was  necessary 
to  simulate  small  disturbances  to  the  projectile  on  the  downward  leg  of 
the  trajectory.  A  small  kick  at  about  50.0  sec  was  imparted  (mathe¬ 
matically)  to  the  round. 

At  about  62  sec  ilight  time  all  supersonic  rounds,  even  those 
for  which  fits  were  not  attempted,  for  one  reason  or  another,  showed  a 
decrease  of  a  few  degrees  in  the  smoothed  level  (not  amplitude)  of  a, 
taking  place  over  one  or  two  precession  periods.  This  could  not  in  any 
way  be  fitted  by  our  model,  and  no  attempt  was  made. 


K,p  had  to  be  used  which  was  grossly  wrong  when  the  behavior 
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1 

1 

I 
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THE  INDIVIDUAL  ROUNDS 


Plates  IA  through  VIE  of  Appendix  B  show  computer  fits  to  both  tra¬ 
jectory  elements  and  yawing  motion  for  rounds  as  follows:  Rd.  11  of 
E-6293,  Rds.  5,  7,  8  and  10  of  E-727»*A  (3  Apr  1968)  and  Rd.  8  of  E-727*»B 
(lU  Kay  1968).  For  each  set  of  plates  the  Roman  numeral  indicates  the 
round.  The  plate  with  the  A  designator  is  a  plot  of  the  measured  tra¬ 
jectory  elements,  and  the  B  overlay  shows  the  computed  elements.  Those 
designated  C  (and  VE)  are  xerographic  copies  ■"'f  the  yaw  sonde  data  as 
plotted  by  the  1108  computer  and  **060  plotter,  while  those  designated  D 
(and  VF  and  VIE)  show  the  computed  fit  to  that  data  and  some  auxiliary 
curves.  Nete  that  the  time  scales  of  Plates  IC  *hrough  VC  have  an 


18 


i 


NWC  TP  4' 


arbitrary  origin.  True  zero  time  is  indicated  by  the  left  hand  vertical 
dashed  line.  Time  origin  is  correct  in  VIC.  Note  also  that  the  sign  of 
(sin  <r)  has  been  reversed  in  plotting. 

The  aerodynamic  parameters  used  in  obtaining  the  computed  curves 
are  listed  in  Table  5  for  the  six  rounds  illustrated  and  for  the  four 
additional  ones  which  were  analyzed  in  detail. 

Comments  follow  on  each  round  for  which  any  yawing  data  was 
obtained. 


SUBSONIC  SHOTS 


E-6293,  Rd,  10:  Standard  Band.  Yawing  record  obtained  from  1  to 
2  and  4  to  53  seconds.  Somewhat  noisy.  Almost  no  change  was  found  in 
the  computed  curve  of  sin  c  on  going  from  the  sharp  shouldered  to 
rounded  shouldered  curve  of  K^B^),  though  the  6-  curve  showed  slightly 
larger,  (but  still  very  small)  persistent  nutations.  From  15  to  25 
seconds  the  nutations  were  almost  good  enough  to  fit  for  phase.  A 
plausible  fit  to  frequencies  of  both  nutation  and  precession  was  ob¬ 
tained  by  lowering  the  ratio  I  /  1^  by  2#  from  that  from  Table  1,  and 


by  1$  from  the  value  in  Table  5.  Note,  however,  that  I  and  1^ 

were  not  actually  known  for  this  round  (nor  Rd.  11).  The  values  of 
Table  1  are  means  for  the  lot.  Steady  precession  of  about  4°. 


E-6295,  Rd.  11:  Standard  Band.  Yaw  record  from  1  to  53  sec. 
Nutations  seem  smaller  than  usual.  This  record  was  noisy  and  was  sub¬ 
jected  to  filtering  before  plotting,  which  may  have  attenuated  the 
nutation  amplitude,  though  the  filter  was  designed  to  avoid  this. 
Steady  precession  of  about  4°.  See  Plate  Set  I. 


E-7274A,  Rd.  4:  Standard  Band:  Yaw  record  from  2  to  54  sec^  From 
4  to  24  sec.  the  amplitude  of  the  precession  is  steady  at  about  4°. 

Then  over  about  5  sec,  it  decreases  to  about  half  size  and  stays  so  till 
impact.  Little,  if  any,  nutations  visible  except  initially. 

B-7274A,  Rd.  5=  Standard  Band:  Yaw  record  noisy,  2  to  6  sec,  then 
good  to  54  sec.  Persistent  small  nutations  visible.  Precession  ampli¬ 
tude  large,  about  5°  from  10  to  21  sec,  then  decreasing  as  in  Rd.  4,  to 
about  5° j  .  See  Plate  Set  II . 

E- 7274a,  Rd,  7:  Flush  Band.  Yaw  record  from  7  to  53  sec.  Nuta¬ 
tions  visible  but  unreliable  for  phase  fitting.  Steady  precession 
amplitude  smaller  than  usual,  about  2?8.  See  Plate  Set  III. 


E-7274A,  Rd.  8:  Flush  Band.  Yaw  record  from  0.5  to  53  sec.  Note 
that  no  trajectory  data  were  available  for  this  round.  It  was  analyzed 
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because  it  was  the  only  nearly  complete  yaw  record  for  a  projectile  of 
its  type  and  speed,  and  because  of  the  unexpectedly  small  precession  of 
Rd.  7.  In  obtaining  the  fit  to  the  sonde  data,  the  trajectory  data  were 
assumed  to  be  the  same  as  for  Rd.  7.  Note  that  the  bursts  of  nutations 
seen  in  the  computed  6^  curve  between  nominal  times  of  6  to  8,  11  to 
28,  and  32  to  3^  seconds  are  machine  artifacts.  In  early  fitting  trials 
for  this  round  it  was  found  that  the  growth  rate  of  the  precessions  was 
unaffected  by  concurrent  variations  of  Kjj  and 

following  approximate  relation  was  obeyed:  0.l6  Hq  +  Tq  ■  constant. 

Here  the  subscripts  indicate  subsonic,  zero  yaw,  standard  air  density 
values.  See  Plate  Set  IV. 

E-727^A,  Rd.  9:  Flush  Band.  Similar  to  Rd.  8.  No  trajectory 
available.  Sonde  record  shorter  than  for  Rd.  8. 

E-727^A,  Rd.  10:  No  Band.  Yaw  record  from  3  to  33  sec.  Nutations 
present,  and  clear  from  3  to  19  sec.  To  fit  both  frequencies  between  3 
and  13  sec,  Ip / 1  t  had  to  be  re(iuced  by  0.6$  from  that  deduced  from 

Table  1,  and  by  1.5$  from  that  of  Table  2.  In  order  to  keep  the 

computed  nutations  large  enough  long  enough,  the  drastic  measure  was 

Kp  to  1.41.  This 

pair  gave  unacceptable  growth  of  the  nutation  for  t  >  13  sec,  but  still 
allowed  the  nutations  to  decay  too  rapidly  between  3  and  6  sec.  On 
reinstating  the  original  (Table  5)  values  of  K^,(0)  (0.094l)  and 

.00),  but  keeping  the  new  value  of  I  /l_,  a  good  fit  to  precession 

P  * 

phase  over  the  whole  record  was  secured,  with  a  only  slightly 

larger.  With  the  original  values  of  I  /  3^  and 

precession  growth  was  unaffected  by  joint  change  of  and  so 

long  as  .17^  H  +  T  =  constant.  See  Plate  Set  V. 

00 

E-727^A,  Rd.  11:  No  Band.  Yaw  record  13  to  33  sec.  Some  indi¬ 
cat  ion-ofr^recessIons-8-to-T3-sec ,  but  in  bad  noise.  Nutations  good 
from  20  to  30  sec.  Fit  to  both  frequencies  required  decreasing  Ip/lj 

by  1.2$  and  by  1.4$.  These  new  values  also  gave  a  good  fit  over 

the  whole  record.  Computed  nutation  amplitude  was  maintained  nearly 
steady  as  in  the  actual  data,  over  the  interval  20  to  30  sec  by  using 
1^(0)  of  .1516  and  a  of  1.43,  but  this  set  of  values  again  gave 

unacceptable  nutation  growth  for  t  >  30  sec. 


it  was  found  that 


taken  of  raising  IC^(O)  to  .1568  and  dropping 


so  long  as  the 
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SUPERSONIC  SHOTS 

Unfortunately,  no  good  yawing  records  of  unmodified  (standard  band) 
projectiles  were  obtained  in  this  series  of  experiments.  On  the  other 
hand,  one  good  record  for  the  standard  round  was  obtained  on  8  November 
1965,  during  test  E-5^15,  which  was  conducted  in  connection  with  RAP 
development.  This  record  for  that  test's  Round  6,  extends  from  15  to  80 
sec  time  of  flight.  It  is  smooth,  with  no  visible  oscillations  until 
about  65  sec,  when  a  small  precession  appears,  then  slowly  grows,  attain¬ 
ing  an  amplitude  of  about  2°  or  a  little  more  by  impact.  In  1965  the 
available  radars  were  unable  to  provide  data  smooth  enough  so  that  good 
velocity  data  could  be  derived. 

E-6295,  Rd,  5:  Standard  Band.  Yaw  record  almost  useless  because 
of  noise,  but  it  contains  a  suggestion  of  precession  Just  before  impact. 

E-6293,  Rd.  7:  Standard  Band.  Yaw  record  stops  at  1*0  sec. 

E-7247B,  Rd.  5;  Standard  Band.  Yaw  record  very  noisy.  Precession 
visible  from  57  sec.  A  shift  in  average  level  of  sin  o-  occurs  at  65 
sec  without  noticeably  changing  the  precessional  amplitude.  (This  type 
behavior  was  found  for  all  rounds  of  E-7274B. ) 

E-7274B,  Rd.  7:  Flush  Band.  Yaw  record  from  4  to  76  sec.  Similar 
to  Rd.  8.  Precession  becomes  visible  about  55  sec.  Shift  at  64  sec. 

Less  trajectory  data  them  for  Rd.  8. 

E- 7274b,  Rd.  8;  Flush  Band.  Yaw  record  from  5  to  76  sec.  Preces¬ 
sion  visible  from  51  sec.  Shift  at  64  sec.  The  1^(0)  value  of  Table  5 

(•1393)  gave  the  fit  of  Plate  VID  when  the  round  was  given  a  small  kick 
at  t  *  47.5  sec.  The  fit  of  Plate  VIE  was  obtained  by  eliminating  the 
kick,  but  increasing  1^(0)  to  O.1857.  1^(0)  of  0.0928  used  with  a 

double- si zed  kick  gave  an  amplitude  of  about  the  right  eventual  size 
but  too  slow  a  growth.  See  Plate  Set  VI. 

E-7274B,  Rd.  10:  No  Band.  Yaw  record  from  15  to  75  sec.  Preces¬ 
sion  visible  from  52  sec.  Shift  at  6i  sec.  Amplitude  at  75  sec  about 
5°.  Got  fit  with  K_(0)  of  0.1552  and  a  small  kick  at  t  =  53.0  sec. 
Also  got  c  good  fit  with  same  K_  without  kick  but  with  large  initial 
yaw.  1 


E-7274B,  Rd.  11:  No  Band.  Yaw  record  from  16  76  sec.  Similar 
to  Rd.  10.  Precession  visible  from  49  sec.  Shift  at  62  sec.  Amplitude 
at  76  sec,  a  little  less  than  4°. 
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DISCUSSION 


The  drag  values  deduced  from  the  trajectory  fits  must  not  he  taken 
too  seriously,  especially  the  very  low  scale  factors  k^  and  the  large 

flat  corrections  found  for  the  two  supersonic  shots.  In  the  first  place, 
the  velocity  determinations  are  among  the  most  imprecise  of  the  data 
obtained.  In  the  second  place,  no  attempt  was  made  to  get  the  best 
possible  drag  values.  So  long  as  fits  V(t)  and  y(t),  good  enough  to 
serve  as  a  basis  for  calculating  the  angular  motion,  were  obtained,  no 
more  was  required.  In  the  third  place,  no  allowance  was  made  for  wind. 
This  neglect  was  probably  particularly  influential  in  affecting  the 
deduced  values  of  drag  for  the  supersonic  shots,  for  in  those  cases 
there  was  a  large  head  wind  at  altitudes  above  9,000  ft,  and  it  is 
primarily,  though  not  only,  the  velocity  behavior  at  high  altitude  that 
required  the  large  flat  correction,  compensated  on  the  upward  trajectory 
leg  by  a  low  k^. 

The  facts  that  in  seven  cases  the  sun  had  to  be  "moved"  slightly, 
and  that  even  so,  the  general  trends  of  the  sin  a  curves  are  not  per¬ 
fect  over  the  whole  flights  are  probably  due  to  two  causes:  slight 
errors  in  sonde  calibration  and  neglect  of  wind. 

Both  Rounds  4  and  5  of  E-727^A  showed  two  levels  of  steady  preces¬ 
sion  amplitude,  but  the  exactly  similar  Rounds  10  and  11  of  E-6293  did 
not.  In  addition,  the  earlier  level  of  Round  5,  E-727^A  was  unusually 
large.  Round  7  of  E-727^A  did  not  show  the  shift  of  amplitude,  but  did 
have  an  unusually  small  precession,  contrary  to  what  wa6  found  for  the 
similar  Rounds  8  and  9  of  the  same  day's  test.  These  discrepancies  may 
represent  real  effects.  Still  queerer  Magnus  behavior  has  been  found 
for  some  rockets.  But  it  should  be  noted  that  it  is  in  Just  these  three 
cases  that  the  largest  adjustments  to  the  sun  angles  were  needed.  For 
Round  4  the  sun  azimuth  had  to  be  adjusted  by  6°.  For  Round  5  the  sun 
elevation  was  changed  by  11°  and  for  Round  7  by  6°.  It  is  considered 
most  likely  that  these  peculiarities  are  all  due  to  shifts  of  the  sonde 
components  between  calibration  and  flight. 

All  the  rounds  of  E-727^B  (supersonic  shots)  showed  a  peculiar 
shift  of  the  level  of  sin  a  about  which  oscillations  were  occurring. 
This  shift  took  place  at  a  little  over  60  seconds  time  of  flight,  when 
the  projectiles  reached  about  11,000  ft  altitude  on  the  descending  leg 
of  the  trajectory.  Now  there  was  a  strong  wind  shear  between  the  alti¬ 
tudes  of  9,000  and  12,000  ft,  the  speed  changing  from  27  ft/sec  at  the 
lower  altitude  to  92  at  the  higher,  and  the  direction  changed  by  nearly 
90  deg.  Furthermore,  the  shift  was  not  observed  on  Round  6  of  E-5^l5. 

It  is  believed  that  the  shift  was  caused  by  the  wind  shear.  This  belief 
was  tested  by  placing  the  problem  on  a  digital  computer.  The  exact  six- 
degree  -of  -freedom  equations  were  used  and  the  parameters  supplied  for 
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Round  10  of  E-7274B.  Two  runs  were  made,  one  with  and  one  without  the 
measured  wind.  The  latter  portion  of  the  "sin  <r"  plot  for  the  case 
with  wind  is  shown  in  Plate  VII.  One  can  see  the  shift  of  mean  value 
taking  place  at  about  64  sec.  (Remember  that  there  is  an  artificial 
reversal  of  sign  between  the  earlier  plates  and  this  one.)  No  such 
shift  was  found  in  the  curve  computed  with  no  wind. 

The  values  of  K^,  the  spin  damping  coefficient,  scatter  a  good 

deal.  In  particular,  there  is  no  explanation  of  why  larger  values  were 
needed  for  the  two  rounds  of  E-6293  than  for  Rounds  4  and  5  of  E-7274A. 
Comparison  among  the  six  values  for  E-7274A  suggests  that  the  value  of 
Ka  may  really  be  smaller  for  the  rounds  with  flush  rotating  bands  than 

for  those  with  either  the  standard,  engraved  band,  or  no  band  at  all. 

The  bottom  of  the  band  groove  of  the  latter  preserved  the  knurling. 

It  appears  certain  that  from  the  Mach  number  of  0.70  to  at  least 
one  of  0.95  the  Magnus  moment  coefficient  changes  sign  extremely  rapidly 
near  a  yaw  angle  of  about  4°,  for  all  three  versions  of  the  projectile. 

In  fact,  the  data  do  not  exclude  an  actual  discontinuity  in  the  Magnus 
moment.  The  level  and  functional  form  of  the  coefficient  at  larger 
angles  was  not  determinable  from  the  data.  Neither  could  the  data  deter¬ 
mine  the  Magnus  moment  at  supersonic  speeds,  but  no  disagreement  was 
found  with  the  conclusion  from  spark  range  experiments  that  supersonically 
K^,  has  a  small  negative  value. 

The  values  in  Table  5  of  K^{0) ,  the  value  of  K_  at  very  small 
yaw  angles  and  subsonic  speeds  must  not  be  considered  as  definitive.  In 
the  first  place,  they  were  determined  by  considering  only  the  behavior 
of  the  precessions,  and  then  under  the  assumption  that  Kp  was  5.0. 

Without  taking  into  account  the  behavior  of  nutation  amplitude,  neither 
K,p  nor  Kjj  can  be  determined  separately.  Besides  there  is  no  reason 
to  believe  that  the  functional  dependence  of  either  coefficient  on  Mach 
number  and  yaw  is  as  simple  as  assumed.  But  some  simple  choice  of  form 
had  to  be  made.  With  each  additional  complexity  admitted,  the  amount 
and  quality  of  the  data  needed  to  fix  parameter  values,  and  the  uncer¬ 
tainty  of  the  latter,  increase.  The  data  on  precessional  growth,  for 
the  two  rounds  (8  and  10  of  E-7274A)  for  which  this  is  most  completely 
visible,  required  that  very  nearly  K_  / A  +  KR  / 6B  =  constant.  Quick 
trials  on  two  other  rounds  showed  that  this  relation  gave  good  results 
there  also.  Now  the  factor  l/6  corresponds  to  a  stability  factor  of 
about  1.8,  and  to  a  ratio  of  nutation  to  precession  frequency  of  about 
5.  These  are  both  values  appropriate  to  the  conditions  during  the  early 
growth  of  the  precessions.  The  ratio  of  frequencies  was  confirmed  by 
actual  counting.  Note  that  for  the  rounds  fired  supersonically  the 
flight  conditions  during  the  time  they  show  noticeable  precessions  were 
such  that  the  factor  l/6  should  be  replaced  with  one  nearer  l/20. 

Now  most,  but  not  all,  of  the  subsonic  shots  showed  small  persistent 
nutations  superimposed  on  the  much  larger  steady  precession.  This  is 
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not  reproducible  with  the  listed  parameters.  On  the  other  band,  the 
nutation  fitting  done  on  Rounds  10  and  11,  especially  the  latter,  where 
the  Mach  number  range  Involved  was  0.73  to  0.79,  indicates  that  a 
1^(0)  of  0.157,  with  a  corresponding  Kp,  is  too  large.  Tue  behavior 

of  the  two  supersonic  shots  led  to  a  choice  of  K^(0)  of  greater  than 


0.13,  with  Kpj  =  5*0.  It  is,  therefore,  believed  that  very  good  esti- 

and  are  0.l4  and  1.8,  respectively,  for 

0.7  <  u  <  0.9.  In  the  case  of  Round  10,  E-7274A,  the  nutations  were 
even  more  slowly  damped  initially,  that  is  for  0.9  <  u  <1.0,  than 
would  be  computed  from  these  estimated  coefficients.  But  the  preces¬ 
sions  are  growing.  It  is  therefore  likely  that  over  the  latter  range 
of  Mach  number  K_  is  even  smaller.  It  may  even  be  negative  over  a 


mates  for  1(^(0) 


small  interval,  though  there  is  no  proof  of  this.  The  negative  values 
of  Kjj  deduced  from  some  spark  range  shots  are  of  uncertain  signifi¬ 
cance,  for  they  result  from  an  entirely  linear  analysis,  and  are,  at 
best,  badly  scattered. 


In  conclusion,  the  5"0/5^  shell,  Mk  4l,  develops  a  definite  limit 
cycle  in  precession,  of  amplitude  about  4°  when  flying  at  speeds  between 
Mach  numbers  of  0.7  and  1.0.  This  shows  up  drastically  on  low  velocity 
shots,  and  more  gently,  in  the  case  of  shots  fired  at  high  velocity,  in 
the  development  of  decided,  though  probably  harmless  precessions  on  the 
downward  legs  of  long  trajectories.  This  behavior  is  undoubtedly  due 
to  the  large  and  very  nonlinear  Magnus  moment  at  subsonic  speeds. 
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TABLE  1.  Physical  Properties  of  Rounds 


Average  length  overall  26.10  in. 


E.  No. 

Rd. 

No. 

m,  lb 

c.g. 

(in  fr. 
tail) 

lb  in2 

11  2 
lb  in2 

Type 

6293 

10 

11 

68.10 

240.9 

2619 

Std 

Std 

7274a 

4 

67.95 

9.34 

240.9 

2549 

Std 

5 

68.71 

9.39 

240.2 

2669 

Std 

7 

67.58 

9.43 

240.2 

2608 

FI  B 

8 

67.95 

9.49 

241.8 

2644 

FI  B 

10 

66.91 

9.56 

235.5 

2651 

No  B 

11 

66.97 

9.58 

234.3 

2619 

No  B 

7274b 

8 

67.95 

9.50 

241.8 

2645 

FI  B 

10 

67.00 

9.58 

234.7 

2629 

No  B 

TABLE  2.  Firing  Conditions 


Test  No. 

Time6 

Q.E. 

MV 

ft/sec 

Sonde 

Record 

Radar 

Coverage 

E-6293 

15  Dec  66 

— 

5 

Std 

1426 

30° 

2466 

poor 

7 

Std 

1445 

30° 

2561 

short 

good 

10 

Std 

1530 

30° 

1230 

good 

good 

11 

Std 

1536 

30° 

1242 

good 

good 

E-7274A 

3  Apr  68 

4 

Std 

1455 

30° 

1227 

good 

good 

5 

Std 

1504 

30° 

1231 

good 

good 

7 

FI  Bd 

1521 

3°° 

1167 

short 

good 

8 

FI  Bd 

1545 

30° 

1187 

good 

none 

9 

FI  Bd 

1555 

3°° 

1180 

fair 

none 

10 

No  Bd 

1607 

30° 

1200 

good 

late 

11 

No  Bd 

l6l4 

30° 

1200 

good 

late 

E- 7274B 

l4  May  68 

5 

Std 

0732 

45° 

24l6 

poor 

good 

6 

Std 

0747 

45° 

2360? 

poor 

useless 

7 

FI  Bd 

0801 

45° 

2456 

good 

fair 

8 

FI  Bd 

081 4 

45° 

2458 

good 

good 

No  Bd 

0832 

45° 

2474 

good 

good 

No  Bd 

0846 

45° 

2434 

good 

short 

a  Times  for  E-6293  and  E-7274A  are  PST;  for  E-7274B  they  are  POT. 
**  The  MV  values  are  instrumental  and  are  not  corrected  back  to 
emergence . 
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TABLE  3*  Meteorological  Data 


5293 

7274a 

7274B 

p.103, 

slug/ft ^  T  deg  C 

p*103, 

slug/ft 3  T  deg  C 

p*103, 

8lUg/ft3 

h.ft.MSL 

T  deg  C 

2185 

2.20  16.6 

2.14  22.1 

2.24 

9*7 

3000 

2.17  13.7 

2.11  17.3 

2.18 

8.9 

4000 

2.10  11.3 

2.06  14.4 

2.11 

7.2 

5000 

2.04  9.8 

2.00  11.5 

2.05 

5*0 

6000 

1.97  8.2 

1.96  8.6 

1.99 

2.8 

7000 

1.91  6.8 

1.90  5*6 

1.94 

0.1 

8000 

1.85  5.3 

1.85  2.6 

1.88 

-2.8 

10000 

1.78 

-8.0 

15000 

1.47 

-9*4 

20000 

1.25 

-19.5 

25000 

1.06 

-29.4 

30000 

0.90 

-42.0 

Winds 

<  3  ft/sec,  surface 

<  22  ft/sec,  surface 

<  7  ft/sec 

,  surface 

<  27  ft/ sec  to 

<  16.8  ft/sec  to 

increasing  to  33  at 

8000  ft 

8000  ft 

11,500  ft  and  to 

90  at  14,000  ft 

Used 

p,  Sur- 

2.21*10° 

2.15-10-5 

2.24*10' 

-3 

face 

l6?0  C 

l4?5  c 

Temp 

0°.  C 

a 

1119  ft  sec 

lll6  ft  sec 

1087.6  ft  sec 

7 

.306*10'^  per  ft 

.258 *10" 4 per  ft 

.316*10' 

’^per  ft 

TABLE  4,  Drag  Function 


Mach  No. 

kdoo 

Mach  No. 

KD00 

0.0 

0.0475 

1.3 

0.1448 

0.7 

0.0475 

1.4 

0.1390 

C.8 

0.0480 

!  1.6 

0.1302 

0.9 

0.0558 

2.0 

0.1170 

1.0 

0.1261 

2.2 

0.1105 

1.1 

1.2 

0.1539 

0.1488 

2.4 

0.1041 
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PLATE  I.  E4293,  Round  1». 

A.  Trajectory  data. 

R.  Computed  trajectory  elements. 


PLATE  I.  (Contd.) 

C.  Yaw  sonde  record. 

D.  Computer  simulation  of  yaw  sonde  record. 


PLATi:  II.  U-7274A,  Round  5. 

A.  Trajectory  data. 

B.  Computed  trajectory  elements. 
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PLATE  H.  (Contd.) 

C.  Yaw  sonde  record. 

D.  Computer  simulation  of  yaw  sonde  record. 
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PLATE  IU.  E-7274A,  Round  7. 

A.  Trajectory  data. 

B.  Computed  trajectory  elements. 
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PLATE  III.  (Comd.) 
t'.  Yaw  sonde  record. 

D.  Computer  simulation  of  yaw  sonde  record. 
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PLATE  IV.  E-7274A,  Round  8. 

A.  Trajectory  dau. 

B,  Computed  trajectory  elements. 
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PLATE  IV.  (Contd.) 

C.  Yaw  sonde  record. 

D.  Computer  simulation  of  yaw  sonde  record. 
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PLATE  V.  E-7274A,  Round  10. 
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PLATE  VI.  (Conld.) 

C.  Yaw  sonde  record. 

D.  Computer  umuUtton  of  yaw  sonde  record;  projectile  disturbed  *1  47, g  sec. 


u 


PLATE  VI.  K 
C.  Yaw  soi 
E.  Comput 
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PLATE  VII.  Computed  Yaw  Sonde  Curve  for  E-7274B,  Round  10,  Showing  Effect  of  Wind  Shear 
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NOMENCLATURE 

Dimensionless  open  decay  coefficient 

Azimuth  of  sun,  measured  from  line  of  fire 

Velocity  of  sound 

Coefficient,  aerodynamic  notation 

Dimensionless  drag  coefficient 

Diameter 

Elevation  of  sun 
Force  component 

Dimensionless  longitudinal  gravity  component 
Dimensionless  transverse  gravity  component 
Acceleration  of  gravity  (scalar) 
Dimensionless  yaw  damping  coefficient 
Polar  moment  of  inertia 
Transverse  moment  of  inertia 
The  unit  imaginary 
coefficient,  ballistic  notation 
Dimensionless  left  coefficient 
Dimensionless  overturning  moment  coefficient 
Moment  component 
MS8S 

Dimensionless  normal  force  coefficient 
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p  *  21m 

1  2 
Q  -  5  pv^ 

q 

r 

S  «  ird2/4 
T 
t 

U  *  V/d 
u 
V 


v 

V 

X 


Spin  rate  (rev/ sec) 

Axial  angular  velocity  component 
Dynamic  pressure  component 
Horizontal  angular  pressure  component 
Third  angular  pressure  component 
Area 

Dimensionless  Magnus  moment  coefficient 
Time 

Dimensionless  velocity 
Axial  velocity  component 
Total  velocity 

Horizontal  velocity  component 
Third  velocity  component 
Horizontal  distance  along  firing  line 


y  Altitude 


a 

y 

8*0-9 
e  *  92  +  ie 
M  *  v/a 
P 
<r 

0  -  02  ♦  i0j 


Particle  trajectory  angle 
in  decrement  of  air  density 
Complex  ya«r  angle 

Complex  trajectory  deviation  angle 
Mach  number 
Air  density 

Orientation  angle  referred  to  sun 
Orientation  deviation  angle 
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